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ABSTRACT 

We report the discovery of HAT-P-25b, a transiting extrasolar planet orbiting the V=13.19 
G5 dwarf star GSC 1788-01237, with a period P = 3.652836 ± 0.000019 d, transit epoch = 
2455176.85173 ± 0.00047 (BJEB), and transit duration 0.1174 ± 0.0017d. The host star has a 
mass of 1.01 ± 0.03 M©, radius of 0.96j;QQ4 i?©, effective temperature 5500 ± 80K, and metalhcity 
[Fe/H] — +0.31 ± 0.08. The planetary companion has a mass of 0.567 ± 0.022 Mj, and radius of 
1.190lQQ5g yielding a mean density of 0.42 ± 0.07gcm~^. Comparing these observations with re- 
cent theoretical models, we find that HAT-P-25b is consistent with a hydrogen-helium dominated gas 
giant planet with negligible core mass and age 3.2 ± 2.3 Gyr. The properties of HAT-P-25b support 
several previously observed correlations for planets in the mass range 0.4 < M < 0.7 A/j, includ- 
ing those of core mass vs. metalhcity, planet radius vs. equilibrium temperature, and orbital period 
vs. planet mass. We also note that HAT-P-25b orbits the faintest star found by HATNet to have a 
transiting planet to date, and is one of only a very few number of planets discovered from the ground 
orbiting a star fainter than V=13.0. 

Subject headings: planetary systems — stars: individual (HAT-P-25, GSC 1788-01237) techniques: 
spectroscopic, photometric 



1. INTRODUCTION 

As more transiting extrasolar planets (TEPs) are dis- 
covered, statistics become significant enough to begin 
looking at bulk properties of exoplanet populations. 
By investigating relationships between stellar, plane- 
tary, a.nd orbital characteri s tics (see, e.g.; Hartman et al.l 

' 2O1O0 . we can 



l20ia iBurrows et all [20071: lEnoch et al.l 



probe the underlying astrophysics that dictate the prop- 
erties we observe. In this paper, we present the dis- 
covery of HAT-P-25b, the 25th TEP found by the 
Hun garian-made Autom ated Telescope Network (HAT- 
Net iHakosiEaD I2OOI survey, orbiting the star also 
known as GSC 1788-01237. HAT-P-25b bolsters the pop- 
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ulation of gas giants in the range 0.4 < M < 0.7 Afj, and 
supports several suggestive correlations between plane- 
tary and stellar parameters. 

As we observe more stars in search of TEPs, we also 
expand the parameter space covered by our discoveries. 
One way in which we do so is by following up candidate 
planets around fainter stars. Wide-field ground-based 
transit surveys like HATNet have extensive photomet- 
ric datasets of faint stars from many different fields, but 
have tended to focus their efforts on the brightest stars 
from each field, as these are easier to follow up. Space- 
based surveys like Kepler forego picking bright stars from 
many different fields in order to take advantage of contin- 
uous monitoring of a single field, and thus end up find- 
ing a large percentage of their candidates around faint 
stars. Faint stars pose new challenges, both for pho- 
tometric detection, and for follow-up to rule out false 
positives and confirm planetary status. With the rich 
datasets already in hand from ground-based surveys and 
the rapidly growing number of faint stars observed from 
space, discovering the nuances of faint star follow-up 
will become increasingly important. To date, HAT-P- 
25 is the faintest star around which HATNet has found 
a planet, and, while not yet close to being prohibitively 
faint, it provides a glimpse at some of the obstacles that 
may become routine in the future. 

HATNet has been one of the main contributors to 
the discovery of TEPs. In operation since 2003, it has 
now covered approximately 14% of the sky, searching 
for TEPs around bright stars (8 < / < 14.0). HAT- 
Net operates six wide-field instruments: four at the Fred 
Lawrence Whipple Observatory (FLWO) in Arizona, and 
two on the roof of the hangar housing the Smithso- 
nian Astrophysical Observatory's Submillimeter Array, 
in Hawaii. 
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The layout of the paper is as foUows. In Section [5] we 
report the detection of the photometric signal and the 
follow-up spectroscopic and photometric observations of 
HAT-P-25. In Section |3] we describe the analysis of the 
data, beginning with the determination of the stellar pa- 
rameters, continuing with a discussion of the methods 
used to rule out non-planetary, false positive scenarios 
which could mimic the photometric and spectroscopic ob- 
servations, and finishing with a description of our global 
modeling of the photometry and radial velocities. Our 
findings are discussed in Section |4l 

2. OBSERVATIONS 

2.1. Photometric detection 

The transits of HAT-P-25b were detected with the 
HAT-10 telescope in Arizona and the HAT-9 telescope 
in Hawaii. The region around GSC 1788-01237, a field 
internally labeled as 259, was observed on a nightly basis 
between 2008 September 15 and 2009 March 16, when- 
ever weather conditions permitted. We gathered 6786 
exposures of 5 minutes at a 5.5 minute cadence. Each im- 
age contained approximately 36,000 stars down to Sloan 
r ~ 14.5. For the brightest stars in the field, we achieved 
a per-image photometric precision of 4 mmag. 
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Fig. 1. — Unbinned (top) and binned (bottom) light curves of 
HAT-P-25 including all 6,786 instrumental Sloan r band 5.5 minute 
cadence measurements obtained with the HAT-9 and HAT-10 tele- 
scopes of HATNet (see the text for details), and folded with the 
period P = 3.6528362 days resulting from the global fit described 
in Section [Sj. The solid line shows the "PIPS" transit model fit to 
the light curve (Section 13.311 . 

The calibration of the HATNet frames was carried 
out using standard photometric procedures. The cal- 
ibrated images were then subjec ted to star d e tectio n 
and astrometry, as described in iPal fc BakosI ()2006D . 
Aperture photometry was performed on each image 
at the stellar centroids derived from the Two Micron 
AU Sky Survey (2MASS; Skrut skie et al.ll2006[) catalog 
and the individual astrometric solutions. The result- 
ing light curves were decorrelated (cleaned of trends) 
using the Externa l Parameter Decorrelation (EPD; see 
iBakos et al.l I2010D technique in "constan t" mode and 
the T rend Filtering Algorithm (TEA; see IKovacs et al.l 
[2005h . The light curves were searched for periodic box- 
shaped signa ls using the Box Least-Squares (BLS; see 
IKovacs et al.l [20021 method. We detected a significant 
signal in the hght curve of GSC 1788-01237 (also known 
as 2MASS 03134450-^2511506; a = 03'^13'"44.51s, 5 = 



+25°11'50.7"; J2000; V=13.19. lDFoege et al.ll2006l) . with 
an apparent depth of ~ 15.1 mmag, and a period of 
P = 3.6528 days (see Figure [J). The drop in bright- 
ness had a first-to-last-contact duration, relative to the 
total period, oi q — 0.0321 ± 0.0005, corresponding to a 
total duration of Pq = 2.817 ± 0.041 hr (see Figure [J). 

2.2. Reconnaissance Spectroscopy 

As is routine in the HATNet project, all candidates 
are subjected to careful scrutiny before investing valu- 
able time on large telescopes. This includes spectro- 
scopic observations at relatively modest facilities to es- 
tablish whether the transit-like feature in the light curve 
of a candidate might be due to astrophysical phenom- 
ena other than a planet transiting a star. Many of these 
false positives are associated with large radial-velocity 
variations in the star (tens of kms^^) that are easily rec- 
ognized. We used two facilities to carry out the recon- 
naissance spectroscop y: the Tilling hast Reflector Echelle 
Spectrograph (TRES: lFuresd[200l on the 1.5 m Tilhng- 
hast Reflector a t FLWO, and the Fibre- fed Echelle Spec- 
trograph fFIES: lFrandsen fc Lindberdll999D on the 2.5 m 
Nord ic Optical Telescope fNOT: iDiupvik fc AndersenI 
l2010f ) at La Palma, Spain. Both of these instruments 
provide high-resolution spectra which, with even mod- 
est signal-to-noise (S/N) ratios, are suitable for deriving 
RVs with moderate precision (< 0.3kms~^) for slowly 
rotating stars. We also use the reconnaissance spectra 
to estimate effective temperatures, surface gravities, and 
projected rotational velocities of the target star. One or 
two well timed spectra can very effectively rule out many 
types of false positives, including F-M binaries, grazing 
eclipsing binaries, or systems of three or more stars, all of 
which can be expected to show velocity variation and/or 
composite spectra at levels detectable by these observa- 
tions. Additional false positive tests are performed with 
other spectroscopic material described in the next sec- 
tion. The reconnaissance observations are summarized 
in Table [TJ Below we provide a brief description of each 
instrument, the data reduction, and the analysis. 

Our flrst spectrum was taken with the medium flber on 
TRES, which has a resolving power of A/AA « 44,000 
and a wavelength coverage of ~3900-8900A. The sec- 
ond spectrum was taken with the medium fiber on FIES, 
which has resolving power of A/AA « 46,000 and a wave- 
length coverage of ^3600-7400 A. The spectra were ex- 
tracted and analyzed ac c ording to the procedures out- 
lined by iBuchhave et ahl ()201 0f). Having monitored the 
lAU radial velocity standard HD182488 over the span 
of about 400 days with both TRES and FIES, we can 
compare the velocity zero-points and stability of each in- 
strument to determine if any correction must be applied 
before combining the datasets. We calculate the mean 
TRES velocity to be -20.821 ± 0.060 km s"^ (rms error) 
when correlations are performed in a single echelle order 
against a synthetic spectrum. For FIES, this number 
is —20. 890 ± 0.046 kms^^. From this, we conclude that 
velocities from TRES and FIES are on the same veloc- 
ity scale (to within the errors), and for the purposes of 
detecting velocity variation due to a stellar companion 
(tens of kms~^), no offset need be applied. 

Based on the reconnaissance spectroscopy, we find that 
the system has rms velocity residuals consistent with 
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no velocity variation within the measurement precision, 
and the observations show no evidence of a composite 
spectrum. From this we conclude that HAT-P-25 has 
no stellar companion. Furthermore, the surface gravity 
found by each instrument is consistent with a dwarf star, 
which reduces the likelihood that the HATNet detection 
is caused by a background blend. 

While we found the reconnaissance velocities to be con- 
sistent with each other to within the measurement pre- 
cision, it is worth noting that those velocities were cal- 
culated by making use of only one order of each spec- 
trum, a small wavelength range (~80 A) surrounding the 
Mg I b triplet. By cross-correlating two spectra against 
each other order-by-order and summing the correlation 
functions, we can utilize the information in a larger wave- 
length range, effectively increasing the S/N and reducing 
the measurement errors. In doing so, we may be able to 
detect a statistically significant velocity variation, which 
would be further evidence for a planetary companion. Of 
course, failure to detect such a variation would not rule 
out a planet; it would simply place an upper limit on any 
companion mass. For the multi-order cross-correlation, 
we used observed templates ~ spectra of the lAU stan- 
dard HD182488 taken on the same night as each HAT-P- 
25 spectrum - to perform the analysis. This allowed us to 
shift each relative velocity to the lAU scale and look for 
orbital motion. We used the wavelength range '^4400- 
6650 A, and the resulting velocities are shown in Table[TJ 
The revised velocities imply Mp = 0.74±0.19 Mj, assum- 
ing a mass of IM© for HAT-P-25, which is consistent 
with the classification from the reconnaissance spectra. 
The error does not take into account systematic uncer- 
tainties due to template mismatch, which could affect 
the relative velocities; or the stability of the instruments, 
which could affect the combination of the relative veloc- 
ities onto the absolute scale. Nonetheless, this result 
lends confidence to the decision to move forward with 
more precious follow-up resources for precise velocities. 

2.3. High resolution, high S/N spectroscopy 

Given the significant transit detection by HATNet, and 
the encouraging results from FIES and TRES that rule 
out obvious false positives, we proceeded with the follow- 
up of this candidate by obtaining high-resolution, high- 
S/N spectra to characterize the RV variations, and to 
refine the determination of the stellar parameters. For 
this we used the HIRES instrument (jVogt et all 119941 ) 
on the Keck I telescope located on Mauna Kea, Hawaii, 
between 2009 December and 2010 February. The width 
of the spectrometer slit was 0'.'86, resulting in a resolving 
power of A/AA « 55,000, with a wavelength coverage of 
-3800-8000 A. 

We obtained 8 exposures through an iodine gas absorp- 
tion cell, which was used to superimpose a dense forest 
of I2 lines on the stellar spec trum and establish an accu- 
rate wavelength fiducial (see lMarcv fc ButleHll993) . An 
additional two exposures were taken without the iodine 
cell, for use as a template in the reductions. Relative 
RVs in the solar system bary c entric frame were derived 
as described by iButler et al.1 ()1996D , incorporating full 
modeling of the spatial and temporal variations of the 
instrumental profile. The RV measurements and their 
uncertainties are listed in Table [2l The period-folded 



data, along with our best fit described below in Section |31 
are displayed in Figure [2l 
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Fig. 2. — Top panel: Keck/HIRES RV measurements for 
HAT-P-25 shown as a function of orbital phase, along with our 
best-fit model (see Table [Sjl- Zero phase corresponds to the time of 
mid-transit. The center-of-mass velocity has been subtracted. Sec- 
ond panel: Velocity O—C residuals from the best fit. The error bars 
include a component from astrophysical jitter (3.5 m s"-*-) added in 
quadrature to the formal errors (see Section l3.3ll . Third panel: 
Bisector spans (BS), with the mean value subtracted. The mea- 
surement from the template spectrum is included (see Section I3.2|l . 
Bottom panel: Relative chromospheric activity index S measured 
from the Keck spectra. Note the different vertical scales of the 
panels. 

In the same figure we show also the relative S index, 
which is a measure of the chromospheric activity of the 
star derived from the fiux in the cores of the Ca II H and 
K lines. This index was compute d following the prescrip- 
tion given by Vaughan, Pre ston fc Wilson! (|l978i) . but 
was not cahbrated to their scale (for details on calcu - 
lation of the relative S index, see lHartman et ahl 120091) . 
We do not detect any significant variation of the relative 
S index correlated with orbital phase; such a correlation 
might have indicated that the RV variations could be due 
to stellar activity, casting doubt on the planetary nature 
of the candidate. There is no sign of emission in the 
cores of the Ca II H and K lines in any of our spectra, 
from which we conclude that the chromospheric activity 
in HAT-P-25 is very low. Furthermore, we estimate from 
the HIRES spectra that on the Mt. Wilson activity scale 
adopted bv lKnutson et al.1 ([2010), the logi?^^ = -4.99 
and Shk = 0.166 (median values). The logR'^j^ value 
depends on o ur estimate oi B — V = .75 using the Toff 
calibration of IValenti fc Fischeri ()2005[ ). 
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TABLE 1 

Reconnaissance Spectroscopy of HAT-P-25. 



Intsrumont 


BJD 


RV^' 




log g* 


V sin j 




(2,454,000+) 


(kms-l) 


(K) 




(kms-i) 


TRES 


1113.87949 


-12.417 ±0.036 


5250 ± 125 


4.0 ±0.25 


6.0 ± 1.0 


FIES 


1115.67533 


-12.608 ±0.036 


5250 ± 125 


4.0 ±0.25 


4.0 ± 1.0 



The velocities reported here are the result of a multi-order cross-correlation against 
the lAU standard HD182488, and have been shifted to the absolute velocity scale 
of Nidcvcr ct al. (2002). The errors are calculated from the rms scatter of individ- 
ual orders in the multi-order analysis. This error estimate does not include possible 
systematic errors due to, e.g., template mismatch or instrument stability. 
^ The stellar parameters T^ff ^ , log , and sin i are the result of correlation of a single 
order against a grid of synthetic spectra with assumed [Fc/H] = 0. Errors quoted are 
half of the grid spacing, but do not take into account systematic errors introduced if 
HAT-P-25 has non-solar metallicity. 



TABLE 2 

Relative radial velocities, bisector spans, and activity 
index measurements of hat-p-25. 



BJD 




fTRV*^ 


BS 


CBS 




(2,454,000+) 


(m s"-') 


(m s -*-) 


(ms-l) 


(ms"-"-) 




1188.91303 


-68.72 


2.58 


113.04 


4.61 


0.39 


1190.94964 


60.74 


2.14 


109.06 


4.41 


0.39 


1190.96797 






107.63 


4.51 


0.37 


1191.92204 






108.46 


7.58 


0.40 


1191.93647 


-47.94 


2.31 


95.72 


4.43 


0.36 


1192.93217 


-49.91 


2.62 


34.02 


7.53 


0.40 


1196.88472 


-6.40 


2.36 


-104.19 


30.34 


0.43 


1197.87232 


69.96 


3.01 


-134.17 


36.62 


0.38 


1250.85004 


-77.40 


3.34 


-188.93 


48.96 


0.33 


1251.84005 


13.68 


2.98 


-140.64 


51.95 


0.32 



Note. — Note that for the iodine-free template exposures 
we do not measure the RV but do measure the BS and S index. 
Such template exposures can be distinguished by the missing RV 
value. 

^ The zero-point of these velocities is arb itrar y. An overall off- 
set 7i.oi fitted to these velocities in Section l3.3l has not been sub- 
tracted. 

^ Internal errors excluding the component of astrophysical jitter 
considered in Section 13.31 

Rela tive chromospheric activity index, not calibrated to the 
scale of lVaughan. Preston fc WilsonI II1978I '). 



2.4. Photometric follow-up observations 

In order to permit a more accurate modeling of the 
light curve, we conducted additional photometric ob- 
servations with the KeplerCam CCD camera on the 
FLWO 1.2 m telescope. We observed two transit events 
of HAT-P-25 on the nights of 2009 Oct 31 and 2010 Jan 
1 (Figure [3]). On 2009 Oct 31, 124 frames were acquired 
with a cadence of 133 seconds (120 seconds of exposure 
time) in the Sloan i band, while on 2010 Jan 1, 190 frames 
were acquired with a cadence of 109 seconds (80 seconds 
of exposure time) in the Sloan i band. 

The reduction of these images, including basic cali- 
bration, astrometry, and aperture photomet ry, was per- 
formed as described by IBakos et al.l (j2010t) . We per- 
formed EPD and TFA to remove trends simultaneously 
with the light curve mod eling (for more details, see 
Section[3l and IBakos et al.i r2010). The final time series 
are shown in the top portion of Figure [3l along with our 
best-fit transit light curve model described below; the 
individual measurements are reported in Table |31 

3. ANALYSIS 



TABLE 3 

HiGH-PRECISION DIFFERENTIAL PHOTOMETRY OF HAT-P-25 



BJD 

(2,400,000+) 



Mag 



0"Mag 



Mag(orig)'' Filter 



55136, 
55136, 
55136, 
55136, 
55136, 
55136, 
55136, 
55136, 
55136, 
55136, 



63175 
63284 
63437 
63608 
63764 
63939 
64095 
64269 
64424 
64600 



0.01928 
0.01476 
0.01455 
0.01541 
0.01963 
0.01613 
0.01992 
0.02231 
0.02184 
0.02230 



0.00295 
0.00144 
0.00145 
0.00144 
0.00142 
0.00142 
0.00139 
0.00145 
0.00133 
0.00136 



11.52150 
11.52080 
11.51890 
11.52220 
11.52500 
11.52420 
11.52700 
11.53230 
11.53140 
11.53070 



Note. — This table is available in a machine-readable form 
in the online journal. A portion is shown here for guidance re- 
garding its form and content. 

^ The out-of-transit level has been subtracted. These magni- 
tudes have been subjected to the EPD and TFA procedures, car- 
ried out simultaneously with the transit fit. 

^ Raw magnitude values without application of the EPD and 
TFA procedures. 
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Fig. 3. — Unbinned instrumental i band transit light curves, 
acquired with KeplerCam at the FLWO 1.2 m telescope on 2009 
Oct 31 and 2010 Jan 1. The lig ht c urves have been EPD- and 
TFA-proeessed, as described in ^ 13.31 Curves after the first are 
displaced vertically for c larit y. Our best fit from the global model- 
ing described in Section 13.31 is shown by the solid lines. Residuals 
from the fits are displayed at the bottom, in the same order as the 
top curves. The error bars represent the photon and background 
shot noise, plus the readout noise. 



3.1. Properties of the parent star 

Fundamental parameters of the host star HAT-P-25 
such as the mass (-M^) and radius (i?*), which are 
needed to infer the planetary properties, depend strongly 
on other stellar quantities that can be derived spec- 
troscopically. For this we have relied on our tem- 
plate spectrum obtained with the Keck/HIRES instru- 
ment, and the an alysis package known as S pectroscopy 
Made Easy (SME: IValenti fc Pisk unov"1 996D, along wit h 
the atomic line database of l^lenti & Fischeii (|2005[ ). 
SME yielded the following initial values and uncertain- 
ties: effective temperature Toff* = 5526 ± 88 K, stel- 
lar surface gravity log g* = 4.55 ± 0.09 (cgs), metallicity 
[Fe/H] = -1-0.30 ± 0.08, and projected rotational velocity 
vsini — 1.4 ± 0.7kms^^. 

In principle the stellar effective temperature and metal- 
licity, along with the stellar surface gravity taken as a 
luminosity indicator, could be used as constraints to in- 
fer the stellar mass and radius by comparison with stellar 
ev olution models . Howe ver, because of reasons described 
in iSozzetti et alj (|2007[ ) , we used the a/i?* normalized 
semi-major axis as luminosity indicator instead of log . g*. 
The a/Ri, quantity is closely related to p*, the mean stel- 
lar density, a nd can be derived d irectly from the transit 
light curves (jSozzetti et al.ll2007D and the RV data (for 
eccentric cases, see Section [531) . This, in turn, allows 
us to improve on the determination of the spectroscopic 
parameters by supplying an indirect constraint on the 
weakly determined spectroscopic value of log g^, that re- 
moves degeneracies. We take this approach here, as de- 
scribed below. The validity of our assumption, namely 
that the adequate physical model describing our data is a 
planetary transit (as opposed to a blend), is shown later 
in Section |3^ 



Our initial values of Teff*, logg*,, and [Fe/H] were 
used to determine auxiliary quantities needed in the 
global modeling of the follow-up photometry and radial 
velocities (specifically, the limb-darkening coefficients). 
This modeling, the details of which are described in 
Section [531 uses a Monte Carlo approach to deliver the 
numerical probability distribution of a/i?* and other fit- 
ted variables. For further details we refer the reader to 
iPall ()2009bD . When combining a/Ri, (used as a proxy for 
luminosity) with assumed Gaussian distributions for Tcff^, 
and [Fe/H] based on the SME determinations, a compar- 
ison with stellar evolution models allows the probability 
distributions of other stellar properties to be inferred, 
including logg*. Here we use the stellar evol ution cal- 
culations from the Yonsei-Yale (YY) series by [Yi et al.l 
([2001). The comparison against the model isochrones 
was carried out for each of 20,000 Monte Carlo trial sets 
(see Section [3.3p . Parameter combinations correspond- 
ing to unphysical locations in the H-R diagram (53% of 
the trials) were ignored, and replaced with another ran- 
domly drawn parameter set. The result for the surface 
gravity, log 5* = 4.48 ± 0.04, is not significantly differ- 
ent from our initial SME analysis. However, we carried 
out a second iteration anyway, in which we adopted this 
value of log and held it fixed in a new SME analysis 
(coupled with a new global modeling of the RV and light 
curves), adjusting only Tcff^, [Fe/H], and vsini. This 
gave Teff* = 5500 ± 80 K, [Fe/H] = -^0.31 ± 0.08, and 
usini = 0.5 ± 0.5kms~^, in which the uncertainties for 
the first two have been increased by a factor of two over 
their formal values to include our estimates of the sys- 
tematic uncertainties. A further iteration did not change 
log g^, significantly, so we adopted the values stated above 
as the final atmospheric properties of the star. They are 
collected in Tabled together with the adopted values for 
the macroturbulent and microturbulent velocities. 

With the adopted spectroscopic parameters the model 
isochrones yield the stellar mass and radius Mi, = l.OlOit 
0.032 Afo and R^ = OMd^o^olrRQ, along with other 
properties listed at the bottom of Table Uj HAT-P-25 is a 
G5 dwarf star with an estimated age of 3.2 ± 2.3 Gyr, ac- 
cording to these models. The inferred location of the star 
in a diagram of a/i?+ versus Teff*, analogous to the classi- 
cal H-R diagram, is shown in Figure 21 The stellar prop- 
erties and their la and 2a co nfidence ellipso ids are dis- 
played against the backdrop of [Yi et ahl (|2001[ ) isochrones 
for the measured metallicity of [Fe/H] = -1-0.31, and a 
range of ages. For comparison, the location implied by 
the initial SME results is also shown (triangle). 

The stellar evolution modeling provides color indices 
that may be compared against the measured values as 
a sanity check. The best available measurements are 
the ne ar-infrared magnitud es from the 2MASS Cata- 
logue HSkrutskie et al.ll2006[) . J2MASS = 11.319 ± 0.018, 
-ff2MASS = 10.908±0.021, and /i:2MAss = 10.815±0.018; 
which we have converted to the photometric system of 
the models (ESQ system) using the transformations by 
[Carpenteil ()2001[ ). The resulting measured color index 
is J — i^T = 0.535 ± 0.011. This value is somewhat 
higher than the value oi J ~ K — 0.46 ± 0.02 pre- 
dicted from the isochrones, and suggests that the star 
may be affected by interstellar reddening. Estimates of 
the total reddening along the line of sight may be ob- 
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Fig. 4. — Model isochrones from IYi et al.l 11200 II ) for the measured 
metallicity of HAT-P-25, [Fe/H]= +0.31, and ages of 0.2, 0.5, 1.0, 
2.0, 3.0, 4.0, 5.0, 6.0, 7.0, and 8.0 Gyr (left to right). The adopted 
values of T^ff^. and a/Ri, axe shown together with their Icr and 2a 
confidence ellipsoids. The initial values of T^Si, and a/R^, from the 
first SME and light curve analyses are represented with a triangle. 

tain ed from the dust maps b y ISchlegel et al.l (llQQSlf^ 
and iBurstein fc Heilesl (flOS^) . The average of these 
two values is E{B -V) = 0.159 ± 0.012. The frac- 
tion of this reddening amount that apphes to HAT-P- 
25 depends on the distance to the object and its Galac- 
tic latitude (see, e.g. [Bonifacio, Monai. & Beers 200(J). 
For the distance we use an initial estimate derived from 
the absolute i^T-band magnitude predicted by the mod- 
els (A/k = 3.44 ± 0.11) and the apparent magnitude 
in the Kg band from 2MASS, which is less affected by 
extinction, and which we convert to the system of the 
isochrones as before. The process is iterated, and re- 
suhs in a final reddening of E{B-V) = 0.131 ± 0.011, a 
de-reddened color of J~K — 0.465±0.014 in good agree- 
ment with the models, and a final distance of 297^13 P'^- 
These values are listed in Table S] 

3.2. Spectral line-bisector analysis 

Our initial spectroscopic analyses discussed in 
Section [2?2] and Section [2731 rule out the most obvious as- 
trophysical false positive scenarios. However, more sub- 
tle phenomena such as blends (contamination by an un- 
resolved eclipsing binary, whether in the background or 
associated with the target) can still mimic both the pho- 
tometric and spectroscopic signatures we see. 

Following Torres et al. (2007), we explored the possi- 
bility that the measured radial velocities are not real, 
but are instead caused by distortions in the spectral 
line profiles due to contamination from a nearby unre- 
solved eclipsing binary. A bisector analysis based on the 
Keck spectra was done as described in §5 of iBakos et al.l 

(Hool). 

The bisector spans (BS) show significant variations, 
however they do not correlate with the stellar RV as 
would be expected if the apparent RV v a,riations were 
due to a blend. Following o ur earlier work ([Kovacs et al.l 
[20T0t iHartman et al.l [20091 ) we investigated the effect of 
contamination from moonlight on the measured bisector 
spans (BS). By adopting the s ame definition of the sky 
contamination factor SCF as in lKovacs et al.l (|2010l ). we 
find a correlation between the SCF and BS (see Figure[5l 

see |http: / /irsa.ipac.caltech.edu/applications/DUST| 



TABLE 4 
Stellar parameters for HAT-P-25 



Parameter Value Source 



Spectroscopic properties 

Teffv, (K) 5500 ± 80 SME'' 

[Fe/H] +0.31 ±0.08 SME 

Dsini (kms-i) . . . 0.5 + 0.5 SME 

iimac (kms-i) . . . . 3.60 SME 

fmic (kms-i) .... 0.85 SME 

log-R^^ -4.99 HIRES'' 

7R.V (kms-i) -12.51 + 0.13 TRES+FIES'^ 

Photometric properties 

V (mag) 13.19 TASS 

y-/c (mag) 0.53 + 0.12 TASS 

J (mag) 11.319 + 0.018 2MASS 

H (mag) 10.908 + 0.021 2MASS 

Ks (mag) 10.815 + 0.018 2MASS 

Derived properties 

M* (Mq) 1.010 + 0.032 YY+a/i?*+SME 'I 

i?* (Rq) 0.959lE5'°^^ YY+a/i?*+SME 

logg* (cgs) 4.48 + 6.04 YY+a/R*+SME 

(Lq) 0.75 + 0.10 YY+a/R^,+SME 

Mv (mag) 5.19 + 0.16 YY+a/iJ*+SME 

Mif (mag,ESO) . . 3.44 + 0.11 YY+a/K^+SME 

Age (Gyr) 3.2 + 2.3 YY+a/K*+SME 

E(_B-y) (mag) . . 0.131 + 0.011 2MASS+YY+a/fl;*+SME'= 

Distance (pc) 2971^3 2MASS+YY+a//?^+SME 



^ SME — "Spectroscopy Made Easy" package for the analysis of high- 
resolution spectra {Va lenti fc Piskunovl [19961) . These parameters rely pri- 
marily on SME, but have a small dependence also on the iterative analysis 
incorporating the isochrone search and global modeling of the data, as de- 
scribed in the text. 

^ See Section [273l 

This velocity is on the absolute scale of INidever et aD i2002r . We have 
increased the error to include to include the rms error of individual orders 
in the multi-order correlation as well as an estimate of possible systematic 
errors introduced by the choice of template spectrum and the stability of the 
instruments. 

YY+a/i?.j,+SME = Based on the YY isochrones (lYTeODIlSoJ) , a/R^ as 
a luminosity indicator, and the SME results. 

^ The total E{B — V) along the line of sight comes from ISchlegel et al.j 
119981) and IBurstein &: Heileg 119821) . The final distance and E{B-V) at 
that distance are calculated iteratively using Mk and Ks (see Section l3.H . 

top panel). When we correct for the relation (Figure [SJ 
bottom panel), the residual BS show no correlation with 
orbital phase. There do appear to be a few discrepant 
BS residuals, although this is the faintest star for which 
we have performed this analysis, and the sky contami- 
nation appears to be especially pernicious in this case, 
both of which contribute uncertainties that are not fully 
characterized. As expected, the outliers correspond to 
the spectra with the lowest S/N and strongest SCF. We 
conclude that the observed BS variations are due to con- 
tamination from scattered moonlight, that the velocity 
variations are real, and that the star is orbited by a close- 
in giant planet. 

3.3. Global modeling of the data 

This section describes the procedure we followed to 
model the HATNet photometry, the follow-up photome- 
try, and the radial velocities simultaneously. Our model 
for the follow-up light curves used analytic formulae 
based on lMandel fc Agol (|2002D for the eclipse of a star 
by a planet, with limb darkening being prescribed by a 
quadratic law. The limb darkening coefficients for the 
Sloan i band were interpolated from the tables bv iClareO 
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Fig. 5. — Top panel: The expected BS variation due to sky 
contamination (SCF) compared to the observed variation. The 
correlation between these values indicates that the observed BS 
variation for this relatively faint target appears to be due to con- 
tamination from scattered moonlight. Bottom panel: The resid- 
ual BS varia tion when corrected for the SCF, compared to orbital 
phase. See ij 13.21 



(|2004') for the spectroscopic parameters of the star as 
determined from the SME analysis (Section 13. The 
transit shape was parametrized by the normahzed plan- 
etary radius p = Rp/Ri,, the square of the impact pa- 
rameter b^, and the reciprocal of the half duration of 
the transit C/R*. We chose these parameters because 
of their simple geometric meanings and the fact that 
these show ne gligible correlations (see lBakos et al.ll201Cll : 
lKippinel[2010l ). Our model for th e HATNet data was the 
simplified "P1P3" version of the iMandel fc Agoll ()2002f ) 
analytic functions (an expansion in terms of L egendre 
polyn omials), for the reasons described in iBakos et al.l 
(|2010D . Following the formalism presented bv lPall (j200£ 
the RVs were fitted with an eccentric Keplerian model 
parametrized by the semi-amplitude K and Lagrangian 
elements k = ecosoj and h = esincj, in which uj is the 
longitude of periastron. 

We assumed that there is a strict periodicity in the 
individual transit times. We assigned the transit num- 
ber Ntr = to the first complete follow-up light curve 
gathered on 2010 Jan 01. The adjustable parameters in 
the fit that determine the ephemeris were chosen to be 
the time of the first transit center observed with HAT- 
Net, Tc, -129, and that of the last transit center observed 
with the FLWO 1.2m telescope, T^fi. We used these 
as opposed to period and reference epoch in order to 
minim ize correlations between parameters (see iPal et ahl 
120081) . Times of mid-transit for intermediate events were 
interpolated using these two epochs and the correspond- 



ing transit number of each event, Ntr- The eight main pa- 
rameters describing the physical model were thus Tc,_ 129, 
T'cOi Rp/R*, b"^, C/^*: ^ = ecosw, and h = esincj. 
Three additional parameters were included that have to 
do with the instrumental configuration. These are the 
HATNet blend factor i^inst, which accounts for possible 
dilution of the transit in the HATNet light curve from 
background stars due to the broad PSF (20" FWHM), 
the HATNet out-of-transit magnitude A/q. HATNet, and 
the relative zero-point 7roi of the Keck RVs. 

We extended our physical model with an instrumental 
model that describes brightness variations caused by sys- 
tematic errors in the measurements. This was done in a 
similar fashion to the analysis presented bv IBakos et al.l 
(1201(1). The HATNet photometry has aheady been EPD- 
and TFA-corrected before the global modeling, so we 
only considered corrections for systematics in the follow- 
up light curves. We chose the "ELTG" method, i.e., 
EPD was performed in "local" mode with EPD coeffi- 
cients defined for each night, and TEA was performed in 
"global" mode using the same set of stars and TEA coef- 
ficients for all nights. The five EPD parameters were the 
hour angle (representing a monotonic trend that changes 
linearly over time) , the square of the hour angle (reflect- 
ing elevation), and the stellar profile parameters (equiv- 
alent to EWHM, elongation, and position angle of the 
image). The functional forms of the above parameters 
contained six coefficients, including the auxiliary out-of- 
transit magnitude of the individual events. The EPD 
parameters were independent for both nights, implying 
12 additional coefficients in the global fit. Eor the global 
TEA analysis we chose 20 template stars that had good 
quality measurements for all nights and on all frames, 
implying an additional 20 parameters in the fit. Thus, 
the total number of fitted parameters was 11 (physical 
model with 3 configuration-related parameters) -I- 12 (lo- 
cal EPD) + 20 (global TEA) = 43, i.e. much smaller than 
the number of data points (322, counting only RV mea- 
surements and follow-up photometry measuremen ts) . 

The joint fit was performed as described in Bako s et al.l 
()2010( ) . We minimized in the space of parameters by 
using a hybrid algorithm, con ibining the down hill sim- 
plex method (AMOEBA; see iPress et al.l I1992D with a 
classical linear least squares algorithm. Uncertainties 
for the parameters were derived applying the Markov 
Chain Monte-Carlo method (MCMC, sec Eord 2006) us- 
ing "Hyperplane-CLLS" chains (Bakos et al. 2010). This 
provided the full a posteriori probability distributions of 
all adjusted variables. The a priori distributions of the 
parameters for these chains were chosen to be Gaussian, 
with eigenvalues and eigenvectors derived from the Eisher 
covariance matrix for the best-fit solution. The Eisher 
covariance matrix was calc ulated ana lytically using the 
partial derivatives given byH (pfMl . 

Following this procedure we obtained the a posteriori 
distributions for all fitted variables, and other quantities 
of interest such as a/i?^. As described in Section [3TT1 
a/i?* was used together with stellar evolution models to 
infer a theoretical value for logQi, that is significantly 
more accurate than the spectroscopic value. The im- 
proved estimate was in turn applied to a second iter- 
ation of the SME analysis, as explained previously, in 
order to obtain better estimates of Teff* and [Fe/H]. The 
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global modeling was then repeated with updated limb- 
darkening coefficients based on those new spectroscopic 
determinations. The resulting geometric parameters per- 
taining to the light curves and velocity curves are listed 
in Table El 

Included in this table is the RV "jitter" . This is a com- 
ponent of assumed astrophysical noise intrinsic to the 
star that we added in quadrature to the internal errors 
for the RVs in order to achieve x^/dof = 1 from the RV 
data for the global fit. Auxiliary parameters not listed in 
the table are: Tc -129 = 2454727.55286 ± 0.00233 (BJD), 
Tc.o 2455198.76874 ±0.00050 (BJD), the blending fac- 
tor Binstr = 0.84 ±0.06, and7rei = 13.5±1.8ms-i. The 
latter quantity represents an arbitrary offset for the Keck 
RVs, and does not correspond to the true center of mass 
velocity of the system, which was listed earlier in Tabled 
(7Rv)- 

The planetary parameters and their uncertainties can 
be derived by combining the a posteriori distributions for 
the stellar, light curve, and RV parameters. In this way 
we find a mass for the planet of Mp — 0.567 ± 0.022 A/j 
and a radius of Rp — LlOOtp osg leading to a mean 
density pp = 0.42 ± 0.07 gcm^'^. These and other plane- 
tary parameters are listed at the bottom of Table |Bl We 
note that the eccentricity of the orbit is not significant: 
k = 0.008±0.012, h = -b.020±0.034 (e = 0.032±0.022, 
w = 271 ± 117°). 

4. DISCUSSION 

Fig. m compares HAT-P-25b to other known TEPs on 
a mass-radius diagram. Nothing unusual about it is 
immediately apparent; its mass is between that of Sat- 
urn and Jupiter and its large radius is typical of other 
short-period gas giant planets. However, this "medi- 
ocrity" gives us a rich sample against which to com- 
pare it. There are currently 19 other known transiting 
exoplanets with 0.4 < M < 0.7 Mj, which is a large 
enough sample to rev eal trends in their physical prop- 
erties. For example, lEnoch et al.l ((2010) note the ap- 
parent inverse relationship between planet radius and 
host star metallicity, which is what would be expected 
if the heavy element content of the planets (core masses) 
is pr op ortional to host star metallicity (jGuillot et. all 
|2006[) . iHartman et al.l ()2010[ ) show a related correla- 
tion, that between stellar metallicit y and t he planetary 
core mass inferred from the Fortncv 'et al.l ((2007) mod- 
els. HAT-P-25b is an important addition to these anal- 
yses, as only X0-2b is more metal-rich in this mass 
range. Fig. |7| shows the location of HAT-P-25b in the 
core mass vs. [Fe/H] relation. We also note that HAT- 
P-25b shows agreement with the previously illustrated 
correl ations of planet rad ius vs. equilibrium temperature 
(e.g.; lEnoch et al.l I20ldl ) and orbital period vs. planet 
mass (which was n oted to be a margi nally statistically 
significant result by iSouthworthI ()2010[ )). 

The characteristics of IIAT-P-25b appear to be mostly 
typical of its peers, although there is one aspect in which 
it is rather unusual for a ground-based transit survey 
- at V = 13.19, it represents the faintest system dis- 
covered by HATNet to date, and is one of only a small 
handful of transiting planets discovered by ground-based 
surveys orbiting stars fainter than V = 13. A histogram 
of host star magnitudes is shown in Fig. |8l The fact that 



TABLE 5 
Orbital and planetary parameters 



Parameter Value 

Light curve parameters 

P (days) 3.652836 ± 0.000019 

Tc (BJD) 2455176.85173 ± 0.00047 

Ti4 (days) 0.1174 ± 0.0017 

T12 = T34 (days) 0.0163 ± 0.0018 

a/_R* 10.461°-^^ 

C/R* 19.75 ±0.19 

flp/fl* 0.1275 ±0.0024 

0.208in?^ 

b = a cos i/R^, 0.456lo o9| 

i (deg) 87.6 ± 0.5 

Limb-darkening coefficients ^ 

ai (linear term) 0.3287 

hi (quadratic term) 0.3039 

RV parameters 

K (ms-i) 74.3 ±2.4 

kjiY" 0.008 ±0.012 

/iRv'^ -0.020 ±0.034 

e 0.032 ±0.022 

u (deg) 271 ± 117 

RV jitter (ms"!) 3.5 

Secondary eclipse parameters 

Ts (BJD) 2455178.698 ± 0.027 

14 0.1138 ±0.0060 

Ts,i2 0.0154 ±0.0018 

Planetary parameters 

Mp (Mj) 0.567 ±0.022 

Rp («J) 1-I90l0;[!»l 

C{Mp,Rp) d 0.15 

Pp (gcm-3) 0.42 ± 0.07 

logsp (cgs) S.OOlHe 

a (AU) 0.0466 ± 0.0005 

Teq (K) 1202 ± 36 

0.044 ±0.003 

(F) (10**crgs-lcm-2) f 4.72 ± 0.58 

^ Tc'- Reference epoch of mid transit that minimizes the correlation 
with the orbital period. It corresponds to Ntr — —6. BJD is 
calculated from UTC. T14: total transit duration, time between 
first to last contact; T12 — T34: ingress/egress time, time between 
first and second, or third and fourth contact. 

^ Values for a quadratic law, adopted from the tabulations by 
IClareij I I2004I) according to the spectroscopic (SME) parameters 
listed in Table |4] 

^ Lagrangian orbital parameters derived from the global modeling, 
and primarily determined by the RV data. 

^ Correlation coefficient between the planetary mass Mp and ra- 
dius Rp. 

The Safronov number is given by — ^ (Vcsc/Vorb)^ — 
(a/Rp){Mp/M^,) (see|H fc Barmanll2007l) . 

^ Incoming flux per unit surface area, averaged over the orbit. 

IIAT-P-25 is faint is significant for a couple of reasons. 
First, this demonstrates that HATNet is capable of de- 
tecting planets orbiting stars with V > 13.0. Second, 
and more importantly, follow-up of fainter host stars will 
become more commonplace as we move forward. It is 
obvious from the early results - see Fig. |5| - and the 
known target lists that space-based surveys will provide 
a glut of planets and planet candidates around fainter 
stars, and ground-based surveys will eventually follow up 
their fainter candidates as well. Systems such as HAT-P- 
25b afford the exoplanet community the opportunity to 
hone our follow-up procedures for optimal performance 
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Fig. 6. — Mass-r adius diagram of known TEPs (small filled squares). HAT-P-25b is shown as a large filled square. Overlaid are 
IFortnev et al.l I I2007I) planetary isochrones interpolated to the solar equivalent semi-major axis of HAT-P-25b for ages of 1.0 Gyr (upper, 
solid lines) and 4 Gyr (lower dashed-dotted lines) and core masses of and 10 (upper and lower lines respectively), as well as isodensity 
lines for 0.4, 0.7, 1.0, 1.33, 5.5 and 11.9gcm~^ (dashed lines). Solar system planets are shown with open triangles. 

on faint stars. For example, HAT-P-25b has already il- 
lustrated that faint stars (and their low S/N) need to 
be treated carefully (if not differently) when it conies to 
characterizing sky contamination. Knowledge like this 
can help inform observing procedures and reduction tech- 
niques in the future. 



10 



Quinn et al. 



CO 



CD 

o 
O 

T3 


_C0 

_o 

CO 

O 



80 

60 
40 
20 

-20 
-40 
-60 
-80 
-100 
-120 



1 

- 


1 1 


1 


1 


1 

• 


1 

- 












• 






• 


• 






- 


• 








- 




• 




• 


o 




• 


• 












• 








- 


• 




• 








• 












• 

1 


1 1 


1 


1 


1 


1 



-0.3 



-0.2 



-0.1 







0.1 
[Fe/H] 



0.2 



0.3 



0.4 



0.5 



Fig. 7. — Calculated planet core mass vs. host star metallicity for planets with 0.4 < M < 0.6 Mj. The location of HAT-P-25b is 
shown by an open circle. The core mass for each planet is determined by linear interpolation within the lFortnev et al] l|2007l ) planet model 
tables for the estimated age, mass, and solar-equivalent semimajor axis of the planet. We adopt an age of 4.0 Gyr or 0.3 Gyr for systems 
with an estimated age greater or less than these limits. Planets with a negative calculated core mass have radii that are too large to be 
accommodated by the models. In this case the core mass is linearly extrapolated from the models, and provides a measure for the degree 
to which the observed radius disagrees with the models. 




V mag 

Fig. 8. — Histogram of host star V magnitudes. Blue: Host stars of transiting planets discovered from the ground. Gray: Host stars of 
transiting planets discovered from space. Red: HAT-P-25. 
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